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Abstract 

We present the results of high-resolution ^^CO (J = 1-0) mapping observations 
with the NRO 45m telescope of the area toward the southern bright arm region of 
M 51, including the galactic center, in order to study the physical conditions of the 
molecular clouds in the arm and the interarm. The obtained map shows the central 
depression of the ^'^CO (J = 1-0) emission, the circumnuclear ring (radius ~ 10" - 
20"), and the spiral arm structure. The arm-to-interarm ratio of the ^^CO (J = 1-0) 
integrated intensity is 2 - 4. We also have found a feature different from that found 
in the ^^CO results. For example, the ^^CO distribution shows a depression in part 
of the spiral arm. The ^^CO/^^CO ratio spatially varies, and shows high values (~ 
20) for the interarm and the central region, but low values (~ 10) for the arm. Their 
values indicate that there is a denser gas in the spiral arm than in the interarm. 

The distribution of the ^^CO shows a better correspondence with that of the Ha 
emission than with the ^^CO in the disk region, except for the central region. We 
found that the ^^CO emission is located on the downstream side of the ^^CO arm, 
namely there is an offset between the ^^CO and the ^^CO as well as the Ha emission. 
This suggests that there is a time delay between the accumulation of gas caused by the 
density wave and dense gas formation, accordingly star formation. This time delay is 
estimated to be ~ 10^ yr based on the assumption of galactic rotation derived by the 
rotation curve and the pattern speed of the M 51 spiral pattern. It is similar to the 
growth timescale of a gravitational instability in the spiral arm of M 51, suggesting 
that the gravitational instability plays an important role for dense gas formation. 
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1. Introduction 



One of the most important aims of research is to develop an understanding of the forma- 
tion mechanism of dense molecular gas in a study of star formation in galaxies, because there is 
a clear correlation between the amount of dense gas and star formation, namely, stars are born 
in dense gas clouds. On the other hand, spiral structures, which are a striking feature of galax- 
ies, are thought to have a correlation with dense gas formation and/or star formation because 
the spiral structures are due to density waves that are expected to accumulate or compress the 
molecular gas. 

Recent observational studies at millimeter and submillimeter wavelengths have revealed 
that there exists molecular gas with high temperature and/or high density in external galaxies, 
especially at galactic centers. These molecular gas properties are greatly different from those of 
our Galaxy, such as Giant Molecular Clouds (GMCs) in the Galactic disk. However, only a few 
attempts to understand the properties of molecular gas have so far been made in a galactic disk 
that is on and between the spiral arms (QiO [cite]cite.key-GB93aGarcia-Burillo et al.(1993a)a; 
030 [cite] cite. key-Kuno95Kuno et al.(1995)), because emission is weaker in the disks than in 
the galactic centers. 

Because ^'^CO (J = 1-0) emission is relatively stronger than other molecular line emis- 
sions, except for ^^CO (J = 1-0), many studies have been made of our Galaxy and external 
galaxies. ^^CO is optically thin and a tracer of denser gas (~ 10^~^ cm~'^) than that traced by 
^^CO(J = 1-0) (~ 10^"'^ cm~^), which is optically thick and surrounds dense molecular gas. 
Also, ^^CO is a good tracer of molecular gas mass through the conversion factor from the inte- 
grated intensity of ^^CO to N^i^ (column density of H2). Consequently, the ratio of ^^C0( J = 1- 
0) /^^C0( J = 1-0) can be used as an indicator of gas density with a low or moderate temperature 
and normal metallicity. The ^^CO(J = 1-0) /^'^CO (J = 1-0) ratios are smaller than ~10 for 
the Galactic plane (|00 [cite]cite.key-Sol79Solomon et al.(1979); 000 [cite]cite.key-Pol88Polk 
et al.(1988)), and for the disk of galaxies, e.g. NGC 891 (000 [cite]cite.key-Saka97Sakamoto 
et al.(1997)). On the other hand, the starburst galaxies have higher ratios (000 [cite] cite. key- 
Aal95Aalto et al.(1995); 000 [cite]cite.key-Kiku98Kikumoto et al. (1998)). We must note that 
the high ratio is seen in the high latitude clouds of our Galaxy which are diffuse cold gas (000 
[cite]cite.key-Pol88Polk et al.(1988)). 

To study the nature of molecular clouds in the central region, the spiral arm and the 
interarm of galaxies and to consider the effect of spiral arms on dense gas formation and star 
formation, we present the results of observations using ^^CO(J = 1-0) and Ha emission toward 
the inner disk of the nearby spiral galaxy M 51. M 51 is a famous grand-design spiral galaxy and 
has been investigated by numerous observations at various wavelengths, e.g., optical, infrared, 
radio, and so on. This galaxy shows two beautiful spiral arms, and kinematic indications of 
density waves also have been found within it. We also note that an oval/bar potential is seen 



2 



Table 1. Adopted parameters of M 51. 



Parameter Value 

Center position (2000.0)" 13^29^52.71'' 

47°11'42.6" 

Morphological type'' Sbc 

Systemic Velocity'^ 463 km s~^ 

Distance'' 9.6 Mpc 

Position angle'^ 170° 

Inclination'^ 20° 

a: 0e]citc.key-TH94Turner and Ho(1994) ([cite]cite.key-TH94Turner and Ho(1994)). 

b: ^]citc.key-dV91de Vaucouleurs ct al.(1991) ([cite] cite. koy-dV91de Vaucouleurs et al.(1991)). 

c: [|[;]cite.kcy-Tully74Tully{1974) {[citc]cite.key-Tully74Tully{1974)). 

d: |a]citc.key-ST74Sandage and Tammann(1974) ([cite] cite. key-ST74Sandage and Tammann(1974)). 

in the central region of M 51 which has an Active Galactic Nucleus (AGN: (^0 [cite] cite. key- 
Tera98Terashima et al.(1998)). The parameters of M 51 are summarized in table |I[ 

2. Observations 

2.1. ^^CO observed with the 45m Telescope at Nobeyama Radio Observatory (NRO) 

We performed observations of ^^CO(J = 1-0) emissions toward the south bright arm 
region of M 51, including the galactic center, between 1994 December and 1995 March with 
the 45-m telescope at Nobeyama Radio Observatory (NRO). The full width at a half-power 
beam (FWHP) was 17" at a rest frequency of ^^CO(J = 1-0) (110.2 GHz), which corresponds 
to 782 pc at a distance from M 51 of 9.6 Mpc. The size of the observed region is 100" x 120" 
(4.6 kpc X 5.5 kpc) and the observing grid has an 11" spacing (see figure 1). The absolute 
pointing accuracy was checked every hour using a SiO maser source, and was found to be 
better than 5" (peak value) throughout the observations. The main beam efficiency was ?7mb = 
0.38 at 110 GHz. We used the 2x2 SIS focal-plane array receiver, which can simultaneously 
observe four positions separated on the sky by 34" each (000 [cite]cite.key-Suna95Sunada et 
al.(1995)). The typical system noise temperatures, including atmospheric effects and antenna 
ohmic losses, were about 400 K. The 2048-channel wide-band acousto-optical spectrometers 
(AOS) were used as receiver backends. The frequency resolution and the channel spacing were 
250 kHz and 125 kHz, respectively, and the velocity resolution was 0.68 km s^^ at 110 GHz. 
The total bandwidth was 250 MHz, corresponding to 680 km s~^. The integration times of 
each observing point were 80 - 100 minutes in most of the observed region, except for some of 
the points. Because those exceptional points were very short (~ 5 min), we did not use them 
to make the total integrated map (see subsection 3.1). 
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2.2. Ha observed with 65cm Telescope at Gunma Astronomical Observatory (GAO) 

2.2.1. Imaging observation 

A Ha image of M 51 was obtained with the 65cm telescope at the Gunma Astronomical 
observatory (GAO) between 2000 December 29 and 2001 January 11. Although the best seeing 
was 2."5, most of the images used for this study had seeing values of between 3" and 3."5, which 
were good enough for comparison with the CO data taken with the NRO-45m telescope. A 
CCD camera assembled by Apogee was attached to the telescope. Installed in this camera was 
a SITe CCD chip of 1024 times 1024 pixels. The pixel size of this chip, 24 /xm, corresponds 
to 0''6 in the sky and so this camera has a field of view (FOV) of 10.'6 square. This FOV was 
wide enough not only to cover the entire galaxy pair (M 51 and NGC 5195), but to provide 
an area for a sky-level evaluation. This camera was cooled down to — 40°C with a water cycle. 
The temperature was stable within 1°C each night. 

Two narrow-band filters were used with a nominal bandpass of 20 A centered at 6584 A, 
and 6620 A, corresponding to the wavelength of Ha and continuum light, respectively. These 
filters were designed to have a boxy-shaped transmittance curve, so that the nebular emission 
could be detected evenly over the galaxy allowing for the rotation of the galaxy. We must note 
that the redshift of M 51 puts Ha at 6573 A. This suggests that the edge of these filters have a 
much higher sensitivity to the receding side than the approaching side. Because the observed 
region of ^^CO corresponds to the receding side, we can say that there was little effect on the 
obtained Ha distribution compared with the ^^CO map. 

We recognize that images from this camera were subject to the so-called memory effect: 
an image has traces of stars that appeared in preceding exposures. Mainly to avoid this effect, 
we drifted the telescope by approximately 30" for each exposure. Because of angular size of 
the bulge of M 51 is much larger than this, we estimated that the memory effect due to the 
bulge was less than the order of 3% of the Ha nebular emission, though the effect would be 
significantly suppressed by the drift of the telescope and the clipping algorithm in the image 
reduction. 

The exposure time was set at 10 min for all exposures. We restricted ourselves to using 6 
images for each filter with good image quality (less atmospheric extinction, better seeing, and, 
more importantly, less contamination from the memory effect). This gives 60 minutes exposure 
for each filter in total. 

2.2.2. Image reduction 

Reduction was carried out with IRAF in the standard manner. Dark frames were com- 
bined. The mean level of the dark current has a monotonic decrease due to the memory effect, 
since it was obtained after illuminating the twilight sky to take flat-field images. This decrease 
introduced some ambiguity in the dark-subtraction. However, this ambiguity in the dark cur- 
rent was almost flat across the FOV, and we assume that local patterns ( i.e. individual H ll 

4 



regions ) of the nebular emission were not badly corrupted. Flat-fielding was made with com- 
bined twilight sky-flat images. The flat-field images are not flat by an order of 5% due to 
inhomogeneous transmittance across the filter. It is not easy at all to remove this inhomogene- 
ity, and although the final image has a residual of a flat-fielding correction, it is of a much lower 
level (~ 10%) than that of the Ha emission of typical H ll regions. Again, this residual is global 
and would not corrupt the small scale patterns of Ha emission discussed in this paper. 

Hampered by ambiguities in the dark-subtraction and the flat-fielding, we made sky- 
subtraction with a constant value evaluated in the area where it is virtually free from light 
from the galaxies (M 51 and NGC 5195) in the broad-band images. Scaled by the mean flux of 
8 stars appearing in images, 6 images of each filter were combined after positional registration. 
The scaling has typically 5 - 10% ambiguity due to the poor photon statistics of the stars in 
the narrow-band images. The positional registration is acceptable within 0''2. Continuum light 
was subtracted from the Ha image, again, scaled by the 8 stellar fluxes which appeared in both 
images. The stellar fluxes are consistent only within 5%. 

The Ha emission in the bright H ll regions has a signal-to-noise ratio of better than 10, 
if the standard deviation in the sky area is adopted for the noise source. Taking into account 
the ambiguity in the subtraction and the scaling, we estimate that the Ha emission may have 
20% ambiguity, but the local patterns (i.e., smaller than the flat-field inhomogeneity) would 
have a much better confidence level. 

Figure |^ shows the resultant Ha image of M 51. No efforts were made to correct the 
extinction in the parent galaxy. This Ha emission was largely consistent with the previous 
studies. A Ha image of M 51 was first carefully examined by |e]cite.key-Ken89bKennicutt 
et al.(1989) ([cite]cite.key-Ken89bKennicutt et al.(1989)), and several studies followed (000 
[cite]cite.key-Rand90Rand and Kulkarni(1990); 000 [cite]cite.key-ThilkOOThilker et al.(2000)). 
Although quantitative a comparison via detailed photometry of individual H ll regions is not 
our aim, the remarkable similarity of the emission of bright and faint H ll regions would support 
the validity of our observations and data reduction. We note that the paucity of Ha emission 
in the center of NGC 5195 is consistent with the spectroscopic study (000 [cite]cite.key-Ho95Ho 
et al.(1995)), giving additional support to our reduction. 

It should also be added that the 3 stars in the observed region have been used to do 
absolute astrometry. The accuracy has been estimated to a few arcsecs, which is enough to 
compare with the ^^CO and ^^CO data. 

3. Results and Discussion 

3.1. Global Distribution of^^CO Emission 
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Fig. 1. Ha narrow-band image of M 51(NGC 5194). Tlie observed points of ^'^CO are indicated by the 
crosses in the figure. The observing grid uses 11" spacing. 
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3.1.1. Spectral line profile 

Figure ^ shows a line-profile map of the ^^CO emission in the observed region. As can 
be seen in this figure, we detected emission in most of the observed region. However, because 
the emission for a part of the south, west, and north-east regions had poor signal-to-noise, we 
didn't use these data to make an integrated intensity map. It should be noted that the emission 
in the central region of M 51 is very weak, though we found strong emissions surrounding the 
central region, namely the bar-end or the beginning of the spiral arms. In particular, the west 
side of the bar-end shows very strong emission. 

Figure ^ also shows the detection of emission from the interarm as well as strong emission 
along the spiral arm. Figure |^ indicates that a typical peak T^b of the emission in the arm 
is ~ 200 mK and that the west bar-end has very strong emission with Tmb higher than 300 
mK. On the other hand, the peak T^b in the interarm is ~ 100 mK. This means that T^b 
in the arm is higher than that in the interarm by a factor of 2-3. We were also able to find 
a larger velocity dispersion in the arm, 30 - 50 km s~^, than that in the interarm, ~ 20 km 
s~^. The larger velocity width in the arm than in the interarm is consistent with previous 
studies of the ^^CO observations, and can be explained by the streaming motion at spiral arms 
caused by density wave shock (000 [cite] cite. key-Kuno97Kuno & Nakai(1997); 000 [cite] cite. key- 
GB93aGarcia-Burillo et al.(1993a)a). 

3.1.2. Spiral structure 

The spatial structures of the ^^CO emission are clearly revealed in the total integrated 
intensity map (figure 3). Figures 0a, b, and c show the total integrated map of ^'^CO emission 
superposed on a color scale of the ^'^CO emission, itself, the ^^CO emission, and the Ha emission, 
respectively. Figure ^a shows a depression in the ^^CO(J = 1-0) emission in the central region 
of M 51. The distribution also has strong concentrations in both bar-ends, especially on the 
west side. Due to the central depression, we find a ring-like structure surrounding the central 
region. This depression was pointed out by previous interferometric observations of M 51 (000 
[cite]cite.key-Matsu98Matsushita et al.(1998)). From figure 0, we have estimated the arm-to- 
interarm ratio using the 13 and 8 line profiles of the arm and the interarm, respectively. We 
add that the profiles which we cannot distinguish between the arm and the interarm due to the 
location at the edge of the arm have been excluded. Their intensities are 4.2 - 12.6 and 1.3 - 
3.6K km s~^ for the arm and the interarm, respectively. The intensities in the arm significantly 
vary while in principle they are increasing with a smaller radius. On the other hand, we can 
not see the large variation in the interarm, though a slightly larger intensity is seen near the 
beginning of the spiral arm. Comparing the profiles in the arm with the near ones in the 
interarm, we estimate the arm-to-interarm ratio of the total integrated intensity of the ^^CO 
(J = 1-0), to be 2 - 4. The high ratio can be seen at the beginning of the arm or the smaller 
radius due to the high intensity seen there. We can also see the arm-to-interarm ratio in the 
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Fig. 2. Line-profile map of tlie ^"^CO emission spectra witli NRO-45m. Tlie velocity resolution of the 
spectra is 5 km s^^. The x and y axes of the map are the right ascension and the declination in units of 
arcsec and (0, 0) in the figure is the galactic center. The vertical and horizontal axes of each spectrum are 
Tmb (K) from to 0.4 K and Vlsr (km s"!) from 300 to 600 km s-\ respectively. 



8 



azimuthal distribution of the ^^CO in the following figure (figure |T^), showing the same ratio 
as that estimated here. This indicates that the ratio virtually depends on the intensities in 
the arm. Although the ratio seems to be similar to the global ratio of the ^^CO (J = 1-0) 
(000 [cite] cite. key-Kuno95Kuno et al.(1995)), the arm-to-interarm ratio of the ^^CO decreases 
with a smaller radius toward the galactic center, and consequently the arm-to-interarm ratio 
of the ^^CO in the observed region is higher than that of the ^^CO. We note that the ratio 
varies spatially in the spiral arm, showing a high ratio at a smaller radius, or rather near the 
beginning of the spiral arm, as mentioned above. 

These morphological structures are globally in agreement with those of the ^^CO emission 
(figure 1^). However, a close look at the distributions of the ^^CO and the ^^CO will reveal the 
difference between them. First, we can see that there is a depression on the ^"^CO spiral arm 
located at (R.A., Decl.) = (13^29™51^ 47°11'00"), unlike in the ^^CO arm. In other words, 
the extension of the spiral arm of the ^^CO from the bar-end is shorter than that of the ^^CO. 
Next, in figure we find an offset between the ^^CO and the ^'^CO arms. This offset seems 
to be obvious, especially at the beginning of the spiral arm and the peak at (R.A., Decl.) = 
(13^29™57'^, 47°10'50"). Taking account of the motion of gas in the frame of galactic rotation on 
the assumption of the existence of a trailing arm, the ^'^CO arm is located on the downstream 
side of the ^^CO arm. On the other hand, in figure |^c which is the ^'^CO map superposed on 
the Ha emission contour, we can see that the depression on the ^^CO arm is located in a region 
where the Ha emission is not detected. Also, the Ha arm is located on the downstream side 
of the ^^CO arm and the ^^CO arm is located along and closer to the Ha arm than the ^^CO 
arm. These observations indicate that the ^'^CO arm shows a closer correspondence to the Ha 
emission than the ^^CO emission in the arm region. Figure § helps to clarify these features. 
This figure is an i? - 6* plot of ^'^CO emission superposed on 12C0 and Ha emissions, where R 
and 6 are the distance from the center and the azimuthal angle for a face-on view, respectively. 
These figures have a range at from R = 30" to 100", covering the spiral arm and not including 
the central region. It makes the existence of the spiral arm clearer. In the figures, we can see the 
depressions both on the ^'^CO and the Ha spiral arms at i? = 50" while there is no prominent 
depression on the ^^CO arm. We can also see that the ^^CO arm shows good agreement with 
the Ha emission at i? = 30" - 50" and 70" - 80" while they show offsets from the ^^CO arm, 
or rather, they are located on the downstream side of the ^^CO arm. On the other hand, at 
R = 60" - 70", the ^^CO has an offset from both the ^^CO and the Ha. This prompts us to 
ask what is the difference between them? The Ha emission at i? = 60" - 70" is located at 
R.A. = 13''29"^52' and Decl. = 47°10'40". It is clear that the distance of this Ha emission 
from the dust lane and the ^^CO arm is more than twice the distance compared with the other 
Ha emission (cf. figure3 in 000 [cite] cite. key- Rand92Rand et al.(1992)). Therefore, we can say 
that this region is not located on the spiral arm, or rather that it has already escaped from the 
arm. This exceptional region is discussed further in subsection 3.4. Therefore, we can say that 
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the distribution of CO is more similar to Ha than is the CO on the spiral arm, except for 
a part of the arm. We conclude that both the ^'^CO and the Ha structures are located on the 
downstream side of the ^^CO structure. 

3.1.3. Radial distribution 

Figure ^ shows the radial distribution of the ^'^CO emission, along with those of the ^^CO 
emission and the Ha emission. This figure indicates the central depression and the ring-like 
structure of the ^^CO emission at a radius of 20" - 30", corresponding to 0.9 - 1.4 kpc. This 
ring-like structure reflects the central depression and the concentrations at the bar-ends, namely 
the beginnings of the spiral arm. This ring-like structure is very different from the distribution 
of the ^^CO, which gradually decreases toward the outside. Also, the ^^CO emission shows a 
rise at a radius of 55", which is coincident with the location of the spiral arm, while the ^^CO 
gradually decreases toward the outside of the galactic disk. These structures of ^^CO are in 
excellent correspondence with those of the Ha emission, except for the galactic center. On 
the other hand, the radial distribution of the ^^CO emission is different from that of the Ha 
emission. 

The previous studies for the kinematics of M 51 reported the positions of the resonances. 
These studies give the positions of the Inner Lindblad Resonance (ILR), which range from 20" 
to 30", and that of the 4/1 resonance of 50" - 60" (|00 [cite]cite.key-Tully74Tully(1974); |00 
[cite]cite.key-GB93bGarcia-Burillo et al.(1993b)b;0OO [cite]cite.key-Aal99Aalto et al.(1999)). In 
the radial distribution shown in figure]^, we find two peaks at radii of 20" and 50", as mentioned 
above. We suggest that the inner ring corresponds to the ILR. Also, the outer peak is located 
on the radii of the 4/1 resonance. The ^^CO shows only a gradual decrease toward an outside, 
and does not show peak structures at the radii of the ILR and the 4/1 resonance. 

3.2. ^"^CO/^^CO Ratio; Molecular Cloud Properties 

In figure ^ the radial distribution of the ratio of the ^^CO-to-^^CO integrated intensities 
(=i?i2/i3) is also shown along with those of the ^^CO, ^^CO, and Ha emissions. This figure 
indicates that the ratio has a range from 5 to 20, and that the central region shows a high ratio 
of > 20. This high ratio in the central region is due to the depression of ^^CO in the central 
region. Although the radial distribution of the -R12/13 gradually decreases with the radius, there 
is a slight increase between the inner ring of the ^'^CO and the outer ring or the spiral arm. This 
suggests the possibility that the -R12/13 is high between the ring and the spiral arm, namely, in 
the interarm. 

The spatial distribution of -R12/13 is clearer in figure ^ and figure |^, which are the super- 
posed color and grey scale images of the ratio on the Ha emission contour. As can be seen in 
these figures, -R12/13 is found to vary spatially. These diagrams indicate that the low -R12/13 of 
typically ~ 10 can be seen at the position of the spiral arm, except for a part of the region, and 
high -R12/13, >20, at the positions of most of the interarm and the central regions. We must also 
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Fig. 3. (a) Total integrated intensity map of ^^CO emission obtained with NRO-45m. The contour 
interval and the lowest contour is 2 K km s~^, corresponding to 2a. The cross indicates the galactic center 
of M 51. (b) superposition of (a) on the ^^CO emission map (contour; ^0 [cite]cite.key-Naka94Nakai et 
al.(1994)). The contour interval and the lowest contour is 10 K km s~^. (c) superposition of (a) on the Ha 
emission map convolved to the 4" beam (contour). The contour interval is 1 in the lowest contour units. 
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Fig. 4. R - 9 plot of the ^'^CO intensity (grey scale) from R = 30" to 100", superposed on those of (a) 
^^CO emission and (b) Ha emission (contours), respectively. The inclination angle has been corrected. 
The contour intervals are the same as figure ^. The units of the horizontal and the vertical axes are arcsec 
and degree, respectively, and 6 is measured counterclockwise from the minor axis. The Ha emission is 
convolved with the 17" beam. 
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Fig. 5. Radial distributions of the ^^CO emission, the ^^CO emission, the intensity ratio R12/13), and 
the Ha emission. Those of the ^^CO and the Ha arc based on the same region of the observed one in 
"CO. 
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Fig. 6. Map of the intensity ratio ( = i?i2/i3 color scale) superposed on the Ha emission (contour). The 
cross indicates the galactic center of M 51. 

note that the high -R12/13 regions have an excellent anti-correlation with the Ha emission except 
for the galactic center. In other words, the regions with high Ru/is such as the interarm show 
little or no Ha emission, while the low -R12/13 is seen in the arm region where the Ha emission 
are distributed. However, it must be noted that the figure shows that there is a part of the arm 
with high -R12/13 and weak Ha emission, which are located at radii from 30" to 40". We add 
that there are local peaks on the spiral arm. However, their contrasts have been significantly 
lower than the arm-interarm contrast. Only the galactic center has both high R12/13 and strong 
Ha emission, which is from the AGN of M 51. 

We summarize the -R12/13 in various regions of our Galaxy and galaxies in table 2. A high 
-R12/13 has been seen in both IR-bright galaxies (000 [cite]cite.key-Aal95Aalto et al.(1995)) and 
starburst galaxies (e.g., M 82: 000 [cite] cite. key-Kiku98Kikumoto et al. (1998)), and in diffuse 
cold clouds, such as high-latitude clouds (000 [cite]cite.key-Knap88Knap & Bowers(1988)). On 
the other hand, the disk of our Galaxy shows a low -R12/13 value (000 [cite] cite. key-Sol79Solomon 
et al.(1979); 000 [cite]cite.key-Pol88Polk et al.(1988)) and the non-active spirals, such as normal 
and poststarburst galaxies, also have intermediate ratios, ~ 10, in the disks (e.g., NGC 891: 000 
[cite]cite.key-Saka97Sakamoto et al.(1997)). -R12/13 in these non-active galaxies is lower than 
those of the starburst galaxies and the high latitude clouds, and higher than that of the disk 
of our Galaxy. Our results show that the -R12/13 values in the central region and the interarm 
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Fig. 7. R - 9 plot of the Rri/is (gi'cy scale) superposed on the Ha emission (contour). The Ha emission 
is convolved with the 17" beam. 

of M 51 are higher than that in the spiral arm and are closer to those of starburst galaxies and 
high-latitude clouds. It is also indicated that in the arm of M 51 the -R12/13 value is slightly 
higher than those of non-active galaxies. This is because -R12/13 measured on the arm of M 51 
may include a part of that in the interarm, which has a high -R12/13 value, due to the narrow 
arm width (~ 10"), which was obtained by previous interferometric observations with a higher 
spatial resolution (e.g., 000 [cite] cite. key- Rand90Rand and Kulkarni(1990)). Consequently, it 
is reasonable to expect that the obtained -R12/13 on the arm may be higher than that in the 
arm only, and is possibly similar to that expected from the normal galaxies. 

In order to understand the physical properties of the molecular gas in the central and 
disk regions of M 51 (and also other galaxies) from the ^^CO and ^^CO data, we used the Large- 
Velocity-Gradient (LVG) calculations (000 [cite] cite. key-Gol74Goldreich and Kwan (1974); 000 
[cite] cite. key-Sco74Scoville & Solomon(1974)) assuming a one-component model. -R12/13 was 
calculated as a function of the H2 number density from 10^ to 10® cm~^, and the kinetic tem- 
perature from 10 to 1000 K. The collision rates for CO molecules are available from |^]cite.key- 
Flo85Flower and Launay(1985) ([cite]cite.key-Flo85Flower and Launay(1985)) (^ 250 K) and 
|e]cite.key-Mck82McKee et al.(1982) ([cite]cite.key-Mck82McKee et al.(1982)) (^ 500 K). We 
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Table 2. -R12/13 in galaxies. 



galaxy 



region 



i?i2/i3 ReniE 



Our Galaxy 



Disk 5.5 - 6.7"''' 

High latitude clouds 40 



Non-active galaxies 
NGC 891 

NGC 1808 
NGC 4945 
NGC 7331 



Center 
Total 



Center 



15.4'' 

6.6^ 

7.8*= 
10 - 19^ 
6.7 - 8.8f Posts 



Active galaxies 
NGC 7479 
M 82 

NGC 4194 
Arp 299 



IR bright galaxies 



Bar 

Central region 

IC694 nucleus 
IC694 disk 
NGC3690 
center 



5 - >30'' 
20 - >70* 

20^' 

60J 

IQi 

5 - 35'= 



M 51 



Center 



Arm 



Interarm 



8-9' 

> 20™ 
9 - 10' 

10™ 

> 15' 
20™ 



a) : 02]cite.key-Sol79Solomon et al.(1979) ([cite] cite. key-Sol79Solomon ct al.(1979)). 

b) : |e]cite.key-Pol88Polk et al.{1988) ([cite]citc.key-Pol88Polk et al.(1988)). 

c) : |e]cite.key-Knap88Knap & Bowers{1988) ([cite]cite.key-Knap88Knap & Bowers(1988)). 

d) : |e]cite.key-Saka97Sakamoto et al.{1997) ([cite]citc.key-Saka97Sakamoto et al.{1997)). 

e) : |e]cite.key-Aalt94Aalto et al.{1994) ([cite]cite.kcy-Aalt94Aalto et al.{1994)). 

f) : |e]cite.key-Berg92Bergman et al.(1992) ([cite]cite.key-Berg92Bergman et al.(1992)). 

g) : P3]cite.key-Isr99Israel and Baas(1999) ([cite] cite. kcy-Isr99Israel and Baas(1999)). 

h) : |e]citc.key-HuttOOHiittemeister ct al.{2000) ([citc]cite.key-HuttOOHiittemeister et al.(2000)). 

i) : ^]cite.key-Kiku98Kikumoto et al. (1998) {[cite] cite. key-Kiku98Kikumoto et al. (1998)). 
j): 0s]cite.key-Aal97Aalto et al.(1997) {[cite] cite. key-Aal97Aalto^et al.(1997)). 

k): |!]cite.key-Aal95Aalto et al.(1995) {[cite] cite. key-Aal95Aalto et al.(1995)). 
1): ^]cite.key-GB93aGarcia-Burillo et al. {1993a) (1993a). 
m): This work 



fixed the molecular abundances to the 'standard' relative abundance as Z(^^CO) = [^^C0]/[H2] 
= 1 X 10-6 (|00 [cite]cite.key-Sol79Solomon et al.(1979)) and [i2CO]/[i3CO] = 50. We also 
fixed the velocity gradient, dv/dr, as a general value of 1.0 km s-^pc"^. The calculated results 
are shown in figure ^ (000 [cite] cite. key-Matsu00Matsushita(2000)). This diagram indicates 
that -R12/13 depends mainly on density when the temperature is not high (<100 K). In the 
Galactic disk, the temperature of typical molecular clouds has a typical value of 10 K. In this 
case, it seems reasonable that -R12/13 is a good tracer for density of molecular clouds, namely 
-R12/13 decreases with the density. Our results show that a low -R12/13 value, smaller than 10, 
is seen in the arm associated with the Ha emission, and the interarm and a part of the arm 
with no associated Ha emission have a high -R12/13 value of ~ 20, as mentioned above. This 
figure indicates that there are diffuse molecular clouds with a density of 2x10^ cm"^ in the 
interarm and the arm with no associated Ha emission, whereas most molecular clouds in the 
arm associated with the Ha emission have a density larger by a factor of two, in the case that 
the same temperatures of molecular clouds in each region, e.g., 10 K, is assumed. Namely, the 
massive star-forming regions traced by the Ha emission are located in the denser clouds with 
a low -R12/13 value. 

On the other hand, the galactic center as well as the interarm region also has a high -R12/13 
value of > 20, which is not consistent with the results obtained by ||e]cite.key-GB93aGarcia- 
Burillo et al.(1993a)(1993a). The reason for this may be the difference in the spatial resolution 
between their (23") and our (17") observations. Since our -R12/13 map shows a very low ratio 
in the bar-end/beginning of the spiral arm, the larger beam can cause contamination with this 
low-ratio region. To confirm this, we convolved our maps to their resolution and measured the 
ratio at the center. The measured ratio is 12 - 13, which is in agreement with their value (8 - 
9) within the errors. As a result, it is reasonable to expect that -R12/13 decreases in their study. 

Figure |^ shows that the high i?i2/i3 values at the central region of M 51 suggests two 
possibilities for the physical condition of molecular gas: one is a low-density (~ 10^ cm"'^) con- 
dition; the other is a high-density (10^^^ cm"^) and high temperature (> 300 K) condition. To 
define the physical conditions of the molecular gas at the center, we compared our results with 
the previously published high resolution observations. The distribution of molecular gas at the 
central region of M 51 is strongly concentrated toward the center; interferometric ^^CO (000 
[cite]cite.key-Saka98Sakamoto et al.(1998); 000 [cite]cite.key-Aal99Aalto et al.(1999)) and HCN 
(000 [cite] cite. key-Kohno96Kohno et al.(1996)) observations clearly show a centrally peaked dis- 
tribution, and suggest disk-like kinematics surrounding the low-luminosity AGN (see also 000 
[cite]cite.key-Sco98Scoville et al.(1998)). On the other hand, the high-resolution ^'^CO observa- 
tions display a lack of ^'^CO emission toward the center (000 [cite]cite.key-Matsu98Matsushita 
et al.(1998)). The HCN/^^CO intensity ratio suggests that the central region is rich in dense 
molecular gas with a density of ~ lO^-^cm"^ (000 [cite]cite.key-Kohno96Kohno et al.(1996)). 
The HCN/^^CO ratio shows a high value (> 3) and suggests that the central region is domi- 
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KINETIC TEMPERATURE [K] 

Fig. 8. Density (vertical axis) and temperature (horizontal axis) dependence of the ^^CO/^'^CO intensity 
ratio (= i?i2/i3; ^0 [cite] cite. key-Matsu00Matsushita(2000)). The solid contours show curves of constant 
i?i2/i3, and the thick solid curve indicates i?i2/i3 of 10.0. 

nated by hot (> 100 K) and dense (~ 10^=^^ cm-^) gas (000 [cite]cite.key-Matsu98Matsushita 
et al.(1998)), and may be correlated with the AGN (^0 [cite]cite.key-Matsu99Matsushita et 
aL(1999)). If our observations are to be consistent with previous ones, and if the density is 
roughly 10^ cm~^, then our value of i?i2/i3 (20) implies a temperature of several hundred K. 
Active galaxies, such as starbursts, also have high i?i2/i3 values as mentioned above, which may 
also be due to the molecular clouds with high density and high temperature. We can say that 
the central region of M 51 has hot, dense clouds as do active galaxies. 

3. 3. Kinematics 

Figure ^ shows the velocity field of the ^^CO emission, which was made by trac- 
ing the peak velocity of the ^^CO line profile. The previous observations of H I,^^CO and 
Ha emissions for the large area provide a position angle of M 51 of 170° (000 [cite] cite. key- 
Tully74Tully(1974)). In the central region of figure ^a, we found a feature difference from 
that expected from the position angle, P.A. = 170°, namely isovelocity contours show a tilt 
from the minor axis of the galaxy. It is reported from near-infrared observations that there 
is an oval potential in the central region of M 51 (000 [cite]cite.key-Pie86Pierce(1986); 000 
[cite]cite.key-TG88Thronson and Greenhouse(1988)). An oval potential causes a distortion of 
the velocity field; this feature in M 51 was found in previous observations (e.g., 000 [cite] cite. key- 
Kuno97Kuno & Nakai(1997)). 
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The isovelocity contours also show S-shape disturbances from pure rotation at the spiral 
arm. This is consistent with previous studies of ^^CO. This could be due to streaming motion 
caused by a density wave, as was already suggested based on previous studies with ^^CO obser- 
vations (000 [cite]cite.key-GB93aGarcia-Burillo et al.(1993a)a; ^0 [cite]cite.key-Kuno97Kuno 
& Nakai(1997)). 

Because the observations give only the velocity component in the line of sight directly, 
we present a Position- Velocity (P-V) diagram along the major axis (P.A. = 170°) to indicate 
the tangential velocity (figure §3). The width of the strip used to make the P-V diagram along 
the major axis is 17". We find that there is a central depression and a concentration at the 
end of the rigid rotation [R ~ 20") in figure ^d. We can also see a velocity shift of 15 km s~^ 
at the position of the spiral arm indicated by arrow in the figure. These correspond to ~ 40 
- 50 km s~^ on the plane of the galaxy. This was pointed out by previous ^^CO observations 
dOO [cite]cite.key-Kuno97Kuno & Nakai(1997); |00 [cite]cite.key-Aal99Aalto et al.(1999)), and 
it was suggested that molecular gas is shocked across the spiral arm. This figure suggests that 
the dense gas which is traced by ^^CO as well as diffuse gas traced by ^^CO is shocked by the 
density wave. 

3.4- Time Delay of Dense Gas Formation and Star Formation 

As mentioned in section 3.2, the ^^CO emission shows good agreement with the Ha 
emission rather than with the ^^CO one, except for in a part of the arms, and they are located 
on the downstream side of the ^^CO. Here, we discuss the effect of the spiral arms, namely 
the density wave on the dense gas and the star-formation mechanism. Before turning to this 
problem, we must consider the exceptional Ha emission region, which shows an offset from the 
^^CO emission. 

Because this region is far from the dust lane and the ^^CO arm, as mentioned in sub- 
section 3.1.2, it seems reasonable to expect that it has escaped from the spiral arm. If the 
star formation is triggered by a shock originating from the density wave, it indicates that this 
region is an "older, or a more evolved star forming region" than those on the spiral arm, since 
the location of the shock is supposed to be the dust lane. We may say that, in the evolved 
star-forming region, the dense gas decreases due to consumption by the star formation, the 
ionization, the photo-dissociation and so on, and the dense gas formation, itself, is suppressed 
because of a decrease in the surrounding gas when it is far from the arm. For this reason, 
it seems reasonable that the above Ha emission region has no associated dense gas which is 
traced by ^^CO. Therefore, when we consider the effect of the density wave on the formation 
mechanism of the dense gas and the following star, it should be an exception. 

In order to investigate the effect of a spiral arm on the nature of molecular gas, we present 
the azimuthal distributions of ^^CO, ^^CO, -R12/13, and Ha at radii of 30" and 45" (figure p!0[) . 
Although it indicates that the ^^CO and ^^CO emissions increase toward the spiral arm from 
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Fig. 9. (a) Velocity field. The contour interval is 10 km s^^ and the thick line indicates 460 km s^^. The 
dashed lines indicate the major and minor axes, (b) P-V diagrams along major axis. The arrow indicates 
the location of the arm. The contour interval is 15 mK and the lowest contour is the 30 mK. 
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the interarm, their peaks are not coincident, namely, the rise of the CO is delayed compared 
with that of the ^^CO. As a result, -R12/13 also shows a decrease across the spiral arm and the 
trend mentioned above is confirmed. The delays are 10-20 deg at radii from 30" to 50" and 
are smaller at larger radii. This suggests that there is a time delay of the formation of dense gas 
traced by ^^CO from the accumulation of molecular gas by the density wave. The small delay 
(in the unit of degree) at a large radius is reasonable, if there is the same time delay between 
the gas accumulation caused by the density wave and the dense gas formation. The difference, 
d, is 400 - 800 pc at a radius of 48", corresponding to 2.2 kpc. We can estimate the time delay 
corresponding to d from the rotation curve. This time delay is given by = rf/Vamn where 
Varm is the vclocity perpendicular to the arm in the frame corotating with the spiral pattern. 
For the adopted pattern speed of Qp = 14 km s~^ kpc~^ (see 000 [cite]cite.key-Kuno95Kuno et 
al.(1995)) and the pitch angle of 20 degrees, we can obtain the value of Varm = 58 km s~^. This 
is the upper limit of the Varm because we assume a circular motion. Then the corresponding 
time delay is r = (0.7 — 1.5) x 10''' yr at 2.2 kpc. 

We must note that the estimated time delay depends on the adopted pattern speed. That 
is, if a larger pattern speed is adopted, Varm becomes small and r becomes large. Actually, 
the pattern speed obtained by previous studies also has a result of 27 km s~^ kpc"^ (000 
[cite]cite.key-GB93bGarcia-Burillo et al.(1993b)b), and may be larger. In the case of larger 
pattern speed, the above estimate of the time delay is the lower limit. We add that r does not 
increase 1.2 times even if the pattern speed doubles. 

It has been known that the Ha emission, indicating massive star formation, is also located 
on the downstream side of the ^^CO emission (000 [cite]cite.key-Vog88Vogel et al.(1988); 000 
[cite] cite. key-Rand92Rand et al.(1992)). As noted above, it shows good correspondence of the 
^^CO with the Ha in the disk of M 51. From these results, we presume that the star formation 
occurred following dense gas formation later than ~ 10^ yr after the accumulation of diffuse 
molecular gas caused by density wave. 

Here, in order to address the origin of the delay of the dense gas and the star formation 
from the compression of the gas due to the density wave, we consider the gravitational instability 
of molecular gas in the spiral arm. Here, we can introduce the Toomre Q parameter to estimate 
the effect of gravitational instability on the mechanism for dense gas formation (e.g. 000 
[cite]cite.key-Ken89Kennicutt(1989)). The previous studies suggest that the Q parameters have 
a correlation with the star-formation activity (000 [cite]cite.key-Kohno01Kohno et al.(2001); 
030 [cite]cite.key-Shio98Shioya et al.(1998); 000 [cite] cite. key- Tosa97Tosaki and Shioya(1997); 
0DO [cite]cite.key-Ken89Kennicutt(1989)). The Q parameter is expressed as Q = Scrit/Sgas, 
where Sgas is the gas surface density and Scrit is critical density for gravitational instabilities, 
and is used as a criteria for local instability in an isothermal thin disk. Scrit is given by 

E„...„^.T4„(j^)(j_^)m,Pc- (1) 
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where and k are the velocity dispersion and the epicychc frequency for the gas disk, respec- 
tively, a is a dimensionless constant (000 [cite]cite.key-Toom64Toomre(1964)). k is expressed 
as 

0.5 



K 




V{r) dV{r) 



(2) 



r dr 

where V{r) is rotational velocity at the radius of r. In the case that we use 10 km s~^ as the 
velocity dispersion and 0.63 as a (^00 [cite]cite.key-Ken89Kennicutt(1989)), we obtain a Sent 
of 84M0pc~^ and 55 Mopc"^ at the radii of 30" and 50", respectively. It should be added that 
we have used the observational rotation curve, which was fitted with the potential model of 
|e]cite.key-Miya75Miyamoto & Nagai(1975) ([cite]cite.key-Miya75Miyamoto & Nagai(1975)), 
as the rotational velocity and the gradient (see 000 [cite] cite. key-Kuno97Kuno & Nakai(1997)). 
Because the surface density of the gas at the spiral arm is ~ 200 Mopc-2 dOO [cite]cite.key- 
Kuno95Kuno et al.(1995)), the Q values are 0.4 and 0.3. This indicates that the gas in the 
arm is gravitationally unstable. Therefore, dense gas formation in the arm may occur via a 
gravitational instability of the gas; consequently the star formation occurs there. 

Next, we compare the above time scale with those of gravitational instabilities. The 
growth time for a gravitational instability can be estimated as follows (e.g., 000 [cite] cite. key- 
Lar87Larson (1987)): 

where cXv and Sgas are the velocity dispersion and the average surface density of gas, respectively. 
We can adopt these values at a radius of 40" as 10 km s~^ and 200 MqPc"^, respectively, from the 
peak of the ^^CO arm (000 [cite] cite. key-Kuno95Kuno et al.(1995)). The estimated timescale 
is ~ 10^ years, which is similar to the above time delay between the ^^CO and the ^'^CO, and 
we can explain the mechanism of dense gas formation by the gravitational instability. 

We must note that an offset between the ^^CO and the HCN has been found in the central 
region of NGC6951 (000 [cite]cite.key-Kohno99Kohno et al.(1999)), which corresponds to ~ 10^ 
yr. HCN is a tracer of dense gas with a density of 10^~^ cm~^, which is higher than ^'^CO. 
This value is also similar to the timescale of gravitational instabilities. It seems reasonable to 
suppose that dense gas formation occurs after ~ 10^~^ yr from the accumulation of gas. As a 
result, star formation also occurs after that. It is likely that gravitational instability plays an 
important role in the mechanism of dense gas and star formation. 

4. Conclusions 

We present the results of ^^CO(J = 1-0) mapping observations toward the southern 
bright arm region of nearby spiral galaxy M 51 carried out with the Nobeyama 45 m telescope. 
The main conclusions are summarized as follows: 
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Fig. 10. Azimuthal distributions of the ^''CO emission, the ^^CO emission and the intensity ratio obtained 
with NRO-45m. and the Ha emission convolved to the 17" beam at radii of 30" and 45". The horizontal 
and vertical axes of each panel are the azimuthal angles (in degree) and their relative values, respectively. 
The gas flows from left to right in this figure. 

1. We find a depression of the ^^CO emission in the central 20" - 30" region, correspond- 
ing to ~1 kpc. There is a ring-hke structure surrounding the central region, which corresponds 
to the bar-end or the beginning of the spiral arm. 

2. The arm structure in the ^^CO is shown in the total integrated intensity map, and the 
arm-to-interarm ratio is 2 - 4. The ratio of the ^^CO in the observed region is higher than that 
of the ^^CO. Although it shows a global correspondence with the ^^CO, a detailed comparison 
shows spatial differences between them, e.g., there is a depression of the ^^CO on the spiral 
arm. 

3. The ^'^CO distribution shows a good correspondence with that of the Ha emission, 
except for in the galactic center, rather than the ^^CO. The ^^CO and the Ha show the de- 
pression at the same position on the spiral arm and are located on the downstream side of the 

4. A velocity shift is detected at the position of the spiral arm, which suggests a streaming 
motion caused by the density wave. 

5. The ^^CO/^'^CO ratio has a range from 4 to >20. The central region and the interarm 
indicate a high value of ~20, while a low value of ~10 is shown in the arm. They indicate 
that there is denser gas in the arm than in the interarm. The densities in the arm and the 
interarm regions were derived to be > 6 x 10^ and 2 x 10^ cm~^, respectively, based on the LVG 
calculation. The high ratio in the central region is due to very hot and dense gas related to 
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AGN in the nucleus of M 51. 

6. The azimuthal distribution at a constant radius shows that the ^^CO emissions are 
located on the downstream side of the ^^CO emission. This indicates that there is a time delay 
from the accumulation of gas caused by shock to the dense gas formation, and the resultant 
star formation. This time delay leads to ~ 10^ yr from the galactic rotation. This timescale 
is similar to the growth time for the gravitational instability. It suggests that gravitational 
instability plays an important role in dense gas formation. 

We are grateful to NRO staff for the operation and the improvement of the 45-m tele- 
scope. We would like to thank Dr. K. Kohno for useful comment and discussion. We appre- 
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of GAO 65-cm. Y.S. thanks the Japan Society for Promotion of Science (JSPS) Research 
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